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Abstract. 

^ . The observed properties of stars and especially the neutrino 

00 ' signal of the supernova 1987A provide an upper limit to the 

■ axion mass, while the age and expansion rate of the universe 

' provide a lower limit. There remains a "window of opportunity" 

. lO^^eV ^ rua ^ 10~^eV, with large uncertainties on either 

' side, where axions could still exist and where they would provide 

^\ . a significant fraction or all of the cosmic dark matter. The 

' current status of this axion window is reviewed. 

o 

^1 1. Introduction 

The interest in axions as a possible dark matter candidate has recently 
I soared thanks to the heroic experimental progress which has led to full- 

5-H I scale searches for galactic axions in Livermore and Kyoto which may 

^ well turn up dark matter axions within the next few years (Fig. 1). There 

is also a noteworthy axion search program in Novosibirsk |^ . 

For such efforts to make sense one needs to understand the "window of 
opportunity" where axions are not excluded by astrophysical and cosmo- 
logical arguments. It is well known that the requirement that stars not lose 
too much energy by axions leads to a lower limit on the Peccei-Quinn scale 
fa which can be translated into an upper limit on the axion mass by virtue 
of the relationship iria = 0.62 eV (10^ GeV//a). It is also well known that 
the non-thermal production in the early universe leads to an upper bound 
on fa (lower bound on rria) lest axions "overdose" the universe. These 
topics are well-covered in a number of reviews Q and books |^ . Some 
recent refinements, however, justify the present update. 

^ E-mail: rafTelt@mppmu.mpg.de 
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Figure 1. Summary of astrophysical and cosmological axion bounds. The 
globular cluster limit assumes an axion-photon coupling corresponding to 
E/N = 8/3 as in GUT models. 



2. Stellar Limits 

2.1. Globular Cluster Stars 

In analogy to neutral pions, axions generically interact with photons ac- 
cording to Ca-y = ffa7E • Ba with = (a/27r/a) {E/N - 1.92 ± 0.08). 
The parameter E/N is a model-dependent fraction of small integers. This 
coupling allows for the axion decay a — > 27 as well as for the Primakoff 
conversion a ^ 7 in the presence of external electromagnetic fields. Be- 
cause charged particles and photons are abundant in the interior of stars 
they are powerful axion sources. 

A novel energy-loss mechanism would accelerate the consumption of 
nuclear fuel in stars and thus shorten their lifetimes. A case where ax- 
ion emission would be efficient and the stellar lifetime is well established 
are low-mass helium- burning stars, so-called horizontal-branch (HB) stars. 
Low-mass red giants have a degenerate helium core {p w lO^gcm"^, 
T Rs 10^ K) so that axion emission is strongly suppressed relative to the 
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cores of HB stars [p « lO'^gcm"^, T « 10* K) whence the number ratio 
of these stars in globular clusters is a sensitive measure for the operation 
of axionic energy losses. The observed number ratios agree with standard 
theoretical expectations to within a few tens of percent implying a limit 
ga-y ^ 0.6x10^1" GeV"^ In GUT axion models where E/N = 8/3 this 
yields ma ^ 0.4 eV (Fig. 1). Our often-quoted "red giant limit" is less 
restrictive because it was based on the statistically less significant number 
of "clump giants" in open clusters 0. 



2.2. Supernova 1987 A 

Being a QCD phenomenon, axions a generically also couple to nucleons 
N by (CAf/2/a)^^7^75V'Ar9'^a where the Cat for protons and neutrons are 
model-dependent numerical coefficients. The most significant limit on the 
axion-nucleon coupling arises from the cooling speed of nascent neutron 
stars as established by the duration of the neutrino signal from the super- 
nova (SN) 1987A. 

Apart from the well-known overall uncertainty of this argument caused 
by the statistical weakness of only 19 observed neutrinos there is a problem 
related to the axion emission rate from the hot and dense nuclear medium. 
In a naive perturbative picture this rate is computed as a bremsstrahlung 
process NN — > NNa with the nucleons interacting by a spin-dependent 
force. The resulting nucleon spin fluctuations act as a source for the emis- 
sion of axions. Like any other bremsstrahlung process, this rate scales 
with the density as p^. However, if one studies axion emission from the 
more general perspective of linear response theory it turns out that the p^ 
scaling is not maintained to arbitrary densities. Rather, the axion emis- 
sion rate does not exceed a certain limit which can be estimated from 
sum rules of the spin-density dynamical structure functions [||, |^ . While 
a detailed calculation of these structure functions is not available, simple 
estimates indicate that in a SN core one is well in the saturation regime 
of the bremsstrahlung process. The previous naive emission rate is thus 
reduced by about an order of magnitude. 

Moreover, a previous often-quoted limit of nia ^ 1 meV was based 
on the generic coupling constants Cn = 0.5 for both protons and neutrons. 
In realistic axion models these couplings can be much smaller. For example, 
in the popular KSVZ model which is representative for the class of hadronic 
axion models the interaction with neutrons nearly vanishes. Altogether the 
SN 1987A limit is reduced to roughly nia ^ 10 meV [^. Inevitably this limit 
involves large uncertainties which are difhcult to quantify. 

Most recently the SN 1987A axion limit was reexamined in a series of 
numerical cooling calculations where the saturation effect was included and 
the coupling constants to neutrons and protons were chosen appropriately 
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for specific models Q. For KSVZ axions the Hmit was found to be about 
8meV, in perfect agreement with the simpfe estimate of Ref. For 
DFSZ axions the hmit varies between about 4 and 12 meV, depending on 
the assumed value of the angle (3 which measures the ratio of two Higgs 
vacuum expectation values. Because this angle is not known one can only 
use the least restrictive number as a conservative limit. In view of the large 
overall uncertainties it is good enough to remember iria ^ 10 meV as the 
current SN 1987A limit for both KSVZ and DFSZ axions. 

If axions interact too strongly they are trapped and contribute to the 
transfer of energy rather than to a direct cooling of the inner SN core. 
Therefore, axions with nia <^ 10 eV cannot be excluded on the basis of the 
duration of the SN 1987A neutrino signal 

However, axions with masses larger than this, i.e., with stronger interac- 
tions, could actually cause a significant contribution to the signal measured 
in the 1MB and Kamiokande II water Cherenkov detectors by the absorp- 
tion on ^^O and the subsequent emission of 7 rays. To avoid too many 
events one can exclude the range 20 eV ^ iria ^ 20keV 13 1. 



2. 3. White Dwarf Cooling 

In certain models axions couple to electrons by (Ce/2/a)'0e7/i75'0e9''a with 
Ce a model-dependent factor of order unity. For most purposes this deriva- 
tive coupling is equivalent to the pseudoscalar structure — i(?aeV'e75V'ea 
with the Yukawa coupling gae = C'eTOe//a which one may parametrize 
by a26 = (5L/4^)/10-'^. 

It was suggested that axion emission with q;26 ~ 0.45 might dominate 
the cooling of white dwarfs such as the ZZ Ceti star G117-B15A for which 
the cooling speed has been established by a direct measurement of the 
decrease of its pulsation period jl4j. Because of this suggestion a new 
bound on 025 was derived by a method similar to the above number counts 
in globular clusters ||l5|. The resulting limit a26 ~ 0.5 is the best direct 
bound on the axion-electron coupling, but it does not quite exclude the 
possibility that axions could play a certain role in white dwarf cooling. 

The most popular example where axions couple to electrons is the DFSZ 
model where Ce = \ cos^ (3 with (3 an arbitrary angle. In this case one may 
use the SN 1987A limits on the axion-nucleon coupling to derive an indirect 
/3-dependent limit on the axion-electron coupling. In this way I infer from 
Ref. [|lll that the largest axion-electron coupling allowed by SN 1987A 
corresponds to a2& ~ 0.08. For typical parameters of old white dwarfs the 
axion luminosity is 0.7q;26 times their photon luminosity j^. Therefore, 
DFSZ axions do not seem to be able to play a significant role in the cooling 
of old white dwarfs. (A novel cooling mechanism may be required if the 
microlensing events of the MACHO collaboration are to be interpreted 
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as halo white dwarfs p6|.) However, axion coohng may be important in 
strongly magnetized white dwarfs where the cyclotron process provides for 
an additional emission channel |17 . 



3. Mass of Dark Matter Axions 

If axions interacted sufficiently strongly {fa ^ 10* GeV) they would have 
come into thermal equilibrium before the QCD phase transition, leading 
to a background sea of invisible axions in analogy to the one expected for 
neutrinos . This parameter range is excluded by the astrophysical argu- 
ments summarized in Fig. 1 which imply that axions interact so weakly that 
they have never come into thermal equilibrium. Still, the well-known mis- 
alignment mechanism will excite coherent oscillations of the axion field . 
In units of the cosmic critical density one finds for the axionic mass density 

flah^ « 1.9 X 4±i(AieV/m,)i "5 q2 ^(q^^ 

where h is the present-day Hubble expansion parameter in units of 
lOOkms"^ Mpc~^. The stated range reflects recognized uncertainties of 
the cosmic conditions at the QCD phase transition and uncertainties in 
the calculation of the temperature-dependent axion mass. The cosmic ax- 
ion density depends on the initial misalignment angle 8^ which could lie 
anywhere between and tt. The function F{Q) encapsules anharmonic cor- 
rections to the axion potential; for an analytic determination see Ref. |^o| . 

If Qi is of order unity, axions with ma = 0(1 /icV) provide roughly 
the cosmic closure density. The equivalent Peccei-Quinn scale fa = 
0(10^^ GeV) is far below the GUT scale so that one may speculate that 
cosmic inflation, if it occurred at all, did not occur after the PQ phase 
transition. If it did not occur at all, or if it did occur before the PQ transi- 
tion with Trchoat > fa, the axion field will start with a different in each 
region which is causally connected at T « /a. Then one has to average 
over all regions to obtain the present-day axion density. 

More importantly, because axions are the Nambu-Goldstone mode of a 
complex Higgs field after the spontaneous breaking of a global U(l) sym- 
metry, cosmic axion strings will form by the Kibble mechanism pT| . The 
motion of these global strings is damped primarily by the emission of ax- 
ions rather than gravitational waves. At the QCD phase transition the 
U(l) symmetry is explicitly broken (axions acquire a mass) causing domain 
walls bounded by strings to form which get sliced up by the interaction with 
strings. The whole string and domain-wall system will quickly decay into 
axions. This complicated sequence of events leads to the production of 
the dominant contribution of cosmic axions. Most of them are produced 
near the QCD transition at T « Aqcd ~ 200 MeV. After they acquire a 
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mass they are nonrelativistic or mildly relativistic so that they are quickly 
redshifted to nonrelativistic velocities. Thus, the string and domain-wall 
produced axions form a cold dark matter component. 

In their recent treatment of axion radiation from global strings, Battye 
and Shellard [2^ found that the dominant source of axion radiation are 
string loops rather than long strings, contrary to what had been assumed 
in the previous works by Davis |2|] and Davis and Shellard . At a given 
cosmic time t the average loop creation size is parametrized as {£) = at 
while the radiation power from loops is P = k/x with the renormalized 
string tension. The loop contribution to the cosmic axion density is [p3| 



where the overall uncertainty has the same source as before. The exact 
values of the parameters a and k are not known; Battye and Shellard 
expect 0.1 < a/K < 1.0. The expression in square brackets is then between 
0.15 and 1.83. 

The proper treatment of axion radiation by global strings has been 
controversial. Sikivie and his collaborators have consistently argued 
that the motion of global strings was overdamped, leading to an axion 
spectrum emitted from strings or loops with a flat frequency spectrum. In 
Battye and Shellard's treatment, wavelengths corresponding to the loop 
size are strongly peaked; the motion is not overdamped. In Sikivie et al.'s 
picture much more of the string-radiated energy goes into kinetic axion 
energy which is redshifted so that ultimately there are fewer axions. The 
cosmic axion density is then of order the misalignment contribution. 

If axions are the dark matter of the universe one may estimate 0.08 < 
flah^ < 0.40. I have assumed that the universe is older than lOGyr, that 
the total matter density is dominated by axions with 0.3 < rJmattcr < Ij and 
that 0.5 < h < 1.0. Including all of the previously discussed uncertainties 
I arrive at a plausible mass range for dark-matter axions of 



as indicated in Fig. 1. Even though Battye and Shellard tend to favor larger 
axion masses for the dark matter, the treatments of both groups imply 
the same lower end for the plausible range. Evidently, the Livermore and 
Kyoto search experiments cover a large fraction of the predicted range, even 
though it would be desirable to extend the search to even larger masses. 

4. Phase Space Distribution 



Qah^ « 88 X 4±i (1 + „ ^ (fieY/ma) 



1.175 



(2) 




6-250 /ieV Sikivie et al., 
6-2000 ^J.eV Battye and Shellard, 



(3) 



The galactic phase-space distribution of axions may well exhibit novel fea- 
tures which are of relevance for the search experiments. Different initial 
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misalignment angles in different causally connected regions lead to density 
fluctuations which are nonlinear from the start. This leads to the forma- 
tion of "axion mini clusters" which may partially survive galaxy formation 
and thus can be found in the Milky Way today For a suitably large 
initial density contrast these clusters can condense into axionic boson stars 
by virtue of higher order axion-axion couplings and may be detectable by 
the femtolensing effect |2^. The direct search experiments need sufhcient 
sensitivity to pick up the diffuse component of the galactic axions which 
are not locked up in mini clusters. 

Another interesting possibility is that axions may have maintained some 
of their initial phase-space distribution, i.e. that they are not fully virialized 
in the galaxy. In this case the axion velocity distribution would exhibit 
very narrow peaks which could enhance the sensitivity of the direct search 
experiments and in any case may be detectable in the laboratory p7[|. 



5. Summary 

After some refinements and corrections, the astrophysical and cosmological 
axion limits seem to have stabilized to what is shown in Fig. 1. There 
remains a "window of opportunity" 10~^ eV ^ ma ^ 10~^ eV where axions 
could still exist. They would then contribute most or all of the dark matter 
of the universe. The ongoing direct search experiments for galactic axions 
have reached a sensitivity where they arc in a position to confirm this bold 
hypothesis or to refute it. 
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